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Abstract 

By using the high- dispersion spectra of 14 bright planet-harboring stars (along 
with 4 reference stars) observed with the new coude echelle spectrograph at Okayama 
Astrophysical Observatory, we investigated the abundances of volatile elements (C, N, 
O, S, Zn; low condensation temperature T c ) in order to examine whether these show 
any significant difference compared to the abundances of other refractory elements 
(Si, Ti, V, Fe, Co, Ni; high T c ) which are known to be generally overabundant in 
those stars with planets, since a T c - dependence is expected if the cause of such a 
metal-richness is due to the accretion of solid planetesimals onto the host star. We 
found, however, that all elements we studied behave themselves quite similarly to 
Fe (i.e., [X/Fe]~ 0) even for the case of volatile elements, which may suggest that 
the enhanced metallicity in those planet-bearing stars is not so much an acquired 
character (by accretion of rocky material) as rather primordial. 
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1. Introduction 

The topic of this paper concerns the chemistry of planet-harboring solar- type stars, which 
have recently been (and are now being) discovered one after another mostly by the Doppler 
technique thanks to the energetic activities of several groups. 

It seems to have been almost established observationally that those stars bearing giant 
planets tend to be mildly metal-rich (i.e., ~+0.2 dex on the average with a considerable scatter) 
compared to the Sun or similar standard stars, though there are some exceptions (see, e.g., 
Gonzalez et al. 2001; Sadakane et al. 2001, Santos et al. 2001; Smith et al. 2001; and 
references therein). What is now in hot debate is, however, the cause of this tendency — why 
does the existence of planets have something to do with the photospheric abundances of host 
stars? 

Here, a promising and often invoked interpretation is the "self-enrichment" hypothesis. 
This attributes the origin of the overabundance of metals to the accretion of large amount of 
H-depleted (i.e., relatively metal-rich) rocky planetesimal materials onto the star, which is more 
or less reasonably understood from the viewpoint of the planet-formation process. 

Following this line, theoretical as well as observational investigations have recently been 
done to clarify whether this mechanism accounts for the observed abundance characteristics 
also in the quantitative sense. Murray and Chaboyer (2001) arrived at a kind of affirmative 
conclusion that adding ~6.5 M e of iron to each star can explain both the mass-metallicity and 
the age-metallicity relations of the stars with planets. Contrary to this, however, Pinsonneault 
et al. (2001) pointed out that an evident T cff -dependence (i.e., more enhanced metallicity 
in earlier-type stars) must be observed if such an accretion mechanism is really responsible, 
which is not observed. Similarly, Santos et al. (2001) concluded by examining the shape of 
the metallicity distribution that the "pollution scenario" due to accretion is unlikely. They 
attributed the higher metallicity of planet-harboring stars to the simple fact that higher gas 
metallicity leads to the high probability of planet formation. 

Meanwhile, another important key to check the "self-enrichment" hypothesis is to study 
the abundance pattern of various elements, especially the difference between the volatile el- 
ements (such as C, N, O) with low condensation temperature T c (i.e., rather unlikely to be 
comprised in rocky materials) and the refractory elements (such as the iron-group elements) 
with high T c (being apt to be condensed into solid particles). Namely, since the elements of 
the former group are expected to be deficient in the accreted materials relative to the latter, 
the abundances of the volatile group would not show such an overabundance as exhibited by 
the refractory group (e.g., Fe). 

Then, is such a trend really observed? Sadakane et al.'s (1999) result for HD 217107, 
that C and O do not appear to share the tendency of marked overabundance shown by the 
other elements, may be suggestive for this prediction. Also, Gonzales and Laws (2000) once 



2 



stated that planet-harboring stars with [Fe/H]> tend to be [C/Fe]< 0, again suggesting that 
C does not keep pace with the enrichment of Fe (though other stars without planets appear 
to show the similar tendency; see also Santos et al. 2000). Later, however, Gonzalez et al. 
(2001) withdrawed this claim by stating "that stars with planets have low [C/Fe] is found to be 
spurious, due to unrecognized systematic differences among published studies" . On the other 
hand, Smith et al. (2001) found that a small subset of planet-harboring stars (i.e., those with 
large planets orbiting quite close) exhibit a trend of increasing [X/H] with an increase in T c , 
which is a sign of the accretion of chemically fractionated solid material. 

Thus, the situation is rather confusing, and a new systematic observational study is 
required to settle this problem. Following this motivation, we decided to investigate, based on 
our own observational data of bright planet-bearing stars collected at Okayama Astrophysical 
Observatory, whether any difference or systematic tendency exists between the abundances of 
volatile and refractory elements, which can be an important touchstone for the adequacy of the 
"self-pollution" scenario. This is the purpose of the present study, in which we will focus on 
the abundances of five volatile elements (C, N, O, S, and Zn) and compare them with those of 
the refractory elements (Si, Ti, V, Fe, Co, and Ni) 

2. Observational Data 

We observed 18 stars (14 planet-harboring stars and 4 reference stars) listed in table 
1 by using the new high dispersion echelle specstrograph (HIDES; cf. Izumiura 1999) at the 
coude focus of the 188 cm reflector of Okayama Astrophysical Observatory in 2000 April, 2000 
October, 2000 December, 2001 February, and 2001 March. 

The slit width was set to 200 fxm (0."76), by which the spectral resolution of R ~ 70000 
was accomplished. By using the single 4Kx2K CCD (13.5 /mi pixel), the wavelength span of 
~ 1200 A could be covered at one time. We have chosen three wavelength regions: green-centered 
region (4900-6100 A; G), red-centered region (5900-7100 A; R), and (near-)IR-centered region 
(7600-8800 A; I). All stars could be observed at both of the wavelength regions G and R, while 
the region I data could not be obtained for six stars (e Eri, 51 Peg, HD 217107, 14 Her, r Cet, 
and (3 Vir). 

The reduction of echelle data (bias subtraction, flat-fielding, scattered-light subtraction, 
spectrum extraction, and wavelength calibration) was performed by using the IRAF 1 software 
package in a standard manner. 

The S/N ratios of the resulting spectra differ from case to case, but are typically of the 
order of 200-300. However, the near-IR spectra turned out to be generally of poorer quality 
because interference fringes could not completely be removed. 

1 IRAF is distributed by the National Optical Astronomy Observatories, which is operated by the Association 
of Universities for Research in Astronomy, Inc. under cooperative agreement with the National Science 
Foundation. 
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Regarding the solar spectrum used as a reference standard, we adopted the very high 
resolution and very high S/N solar flux spectrum published by Kurucz et al. (1984). 

3. Abundance Analysis 

3.1. Model Atmospheres 

With regard to the atmospheric parameters (T eS , logg, [Fe/H], and £) of the program 
stars necessary for constructing model atmospheres and determining abundances, various pub- 
lished studies were consulted. We compared the published values from several sources (e.g., 
Gonzalez 1998; Gonzalez et al. 2001; Gonzalez, Law 2001; Fuhrmann 1998) for a same star 
when available, and found that they are generally consistent with each other and that the 
differences are in most cases within AT c g ~100 K and A logg ~0.1. The adopted parameter 
values (with the references) are presented in table 1. 

We interpolated Kurucz's (1993) grid of ATLAS9 model atmospheres in terms of T e g, 
logg, and [Fe/H] to generate the atmospheric model for each star. 

3.2. Determination of the Abundances 

3.2.1. Equivalent width analysis 

Regarding the target volatile elements (C, N, O, S, and Zn), as far as the lines of interest 
are neither blended nor too weak, we determined the abundances from the equivalent widths by 
using the WIDTH9 program, a companion program to the ATLAS9 model atmosphere program 
written by Dr. R. L. Kurucz (Kurucz 1993). As for the data of spectral lines (wavelengths, 
gf values, and damping parameters), we exclusively invoked the compilation of Kurucz and 
Bell (1995). The gf values of important lines relevant to our analyses are presented in table 2, 
though the choice of absolute oscillator strengths is not essential because the analysis of this 
study is differential relative to the Sun. The observed equivalent widths (which were measured 
by the method of Gaussian fitting or direct integration depending on cases) and the resulting 
abundances are given in tables 3 and 4, respectively. 

Note that the S abundance from the S I 6052 line was determined by correctly taking into 
account the blending two components, while the equivalent width was measured as if regarding 
it as a single line. Also, the non-LTE correction (amounting to ~0. 1-0.4 dex) was taken into 
consideration for deriving the O abundance from the strong O I 7771, 7774, and 7775 triplet 
lines according to the calculation of Takeda (2001), as given in the 7th column of table 4. 

3.2.2. Spectral synthesis analysis 

Meanwhile, refractory element abundances (Si, Ti, V, Fe, Co, and Ni) were determined 
by the spectrum synthetic technique while applying the automatic fitting algorithm (Takeda 
1995), in which the abundances of these 6 elements (along with the macrobroadening parameter 
and the radial velocity shift) were simultaneously varied to find the best-fit solution such that 
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minimizing the O — C deviation. This yields accurate and reliable results than the case of 
conventional analysis. After several trials, the 6080-6089 A region was found to be adequate 
(i.e., a satisfactory fit is accomplished with our gf values given in table 2), which we decided 
to adopt. Note that, apart from those given in table 2, all lines listed in Kurucz and Bell's 
(1995) compilation were included as background opacities. The best-fit theoretically synthesized 
spectra along with the observed spectra are shown in figures la and b. Such established 
abundances are given in table 4. 

Similarly, the C and O abundances from [C i] 8727 and O I 6158 lines were obtained by 
this technique, since these line features are delicately weak and more or less blended with other 
species. As a result of the analyses of these two line features, the abundances of Si and Fe were 
simultaneously obtained as by-products. Note that spectral portions, where the line data are 
so insufficient that the observed spectrum could not be reproduced (e.g., around 6157.3 A for 
the case of lower temperature stars), were masked (i.e., excluded from judging the goodness 
of fit) in the process of automatic fitting. Figures 2 (8727 A region analysys) and 3 (6157 A 
region analysis) show the appearance of the fit, and the resulting abundances are presented in 
table 4. 

3.2.3. Sensitivity to parameter changes 

In order to estimate the uncertainties involved with the ambiguities in the atmospheric 
parameters, we selected v And and 70 Vir as the respresentatives of higher T c fj and lower T e g 
stars and evaluated the abundance variations in response to changing T e s (by +150 K), logg 
(by +0.3), and £ (by +0.5 km s" 1 ). The results are given in table 5, from which we may state 
that errors caused by uncertainties in these parameters are not very significant considering the 
reasonable consistency of the published values (cf. subsection 3.1). 

3.2.4- Final results 

In table 6, we give the final results of the abundances which will be discussed in section 

4. 

As for the volatile elements, we adopted a mean value when two or more abundances 
from different line indicators are available. In this averaging process, a half weight was assigned 
to the C and O abundance derived from the forbidden lines ([C i] 8727 and [O i] 6300) in view 
of their comparatively less reliability due to the suspected blending effect, i.e., with unknown 
species for [C i] (Gonzalez et al. 2001) and with the Ni I line for [O i] (Allende Prieto et al. 
2001). 

Meanwhile, regarding the refractory elements (Si, Ti, V, Fe, Co, and Ni), we adopted 
the abundances from the 6080-6089 A fitting analysis, since we consider the results for these 6 
elements from this 6080-6089 A region is sufficiently accurate. 

The T c values presented in table 6 were taken from Field (1974) for CNO, and from 
Wasson (1985) for the remaining elements. 
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4. Discussion 



4-1. Fe abundance 

The resulting Fe abundances of the program stars (table 6) are compared with the 
literaure values (table 1) in figure 4, where we can see that both are in fairly reasonable 
agreement. The mean for those 14 stars with planets is ([Fe/H]) = +0.12(±0.20), which confirms 
the mild enrichment of metallicity as has been already reported (cf. section 1), though a rather 
large diversity is observed (e.g., p CrB is evidently metal- deficient, though its companion may 
actually not be a planet as remarked in the footnote to table 1). 

4-2. Refractory Elements 

Figures 5a-e show the comparison of the abundances of five refractory elements (Si, Ti, 
V, Co, and Ni) with the Fe abundance. We may state from this figure that the abundances of 
those elements scale with that of Fe, which is an expected tendency because these belong to 
the same group in terms of T c (1300-1500 K). 

The reason for the appreciable systematic difference of [Co/H] by 0.1-0.2 dex relative 
to [Fe/H] (cf. figure 5d) may be related to the weakness of this line (especially for higher T e g- 
stars; cf. figure la), by which the reliability of the results may be more or less lowered. 

4-3. Volatile Elements 

The results of [C/H], [N/H], [O/H], [S/H], and [Zn/H] are compared with [Fe/H] in fig- 
ures 6a-e, respectively. These figures suggest that the abundances of all these five volatile ele- 
ments scale with that of Fe. 2 Especially, the trend of [C/Fe] < for metal-rich planet-harboring 
stars, which was once reported before (cf. section 1), can not be confirmed. Hence, our results 
suggest that the tendency predicted from the "self-enrichment" hypothesis, [X/H] < [Fe/H] 
for metal-enriched planet-harboring stars ([Fe/H] > 0), is not observed at all. (Even a weak 
inversed trend, [X/H] > [Fe/H], appears to exist such as for N or S.) 

If fractionated solid materials have something to do with any chemical peculiarity in 
photospheric abundances, some kind of anti-correlation between the abundances of voilatile 
and refractory elements is to be expected. Let us consider, for example, the case of A-type 
stars (including A Boo stars). In this class of stars, [C/H] is known to be clearly anti-correlated 
with [Si/H], which is suspected as due to the differential accretion of gas and dust in late stages 
of star formation (see, e.g., Holweger, Stiirenburg 1993). On the contrary, in the present case, 
C and Si behave themselves quite similarly with each other, as shown in figure 7 (compare this 
figure with their figure 2). 

2 The trend of appreciable deviation for r Cet (observed also in figure 5) is (at least partly) attributed to 
the consequence of the Galactic chemical evolution for a-process elements (e.g., O, S, Si, and Ti), for which 
metal-deficient stars ([Fe/H] < 0) are known to show [a/Fe] > 0. However, the reason why this tendency 
more or less appears in other (non-a) elements is not clear. 
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We also investigated if there is any systematic trend in [X/H] with respect to T c , such 
that reported by Smith et al. (2001) for a subset of planet-harboring stars. We calculated 
the slope (a) of the linear-regression line ([X/H]=a • T c + b) determined from the [X/H] values 
of the 11 elements (or 10 elements when N is not available). The results are given in 13th 
and 14th columns of table 6. By inspecting these a values, the signs of which are randomly 
positive or negative without showing any systematic dependence with the orbital period/radius 
of the planets, we may state that the tendency of positive a expected from the self-enrichment 
hypothesis (cf. section 1) can not be observed. 

Consequently, we consider that the "pollution" scenario is rather unlikely for the expla- 
nation of the metal-richness of planet-harboring stars, which may lend support for an alternative 
interpretation that the enhanced metallicity of stars with planets is nothing but a "primordial" 
characteristics. 

5. Conclusion 

We determined the abundances of five volatile elements (C, N, O, S, and Zn; low T c ) and 
six refractory elements (Si, Ti, V, Fe, Co, and Ni; high T c ) for 14 bright planet-bearing stars 
along with 4 reference stars (+ Sun), based on the spectra obtained by using the high-dispesion 
echelle spectrograph at Okayama Astrophysical Observatory, in order to examine whether any 
systematic difference exists between these two groups of elements. 

Our results indicate that all of the studied elements behave themselves quite similarly 
with each other (i.e., [X/Fe]~ 0) irrespective of whether the element is volatile (apt to be 
remained in gas phase) or refractory (apt to be included in solid materials). In other words, 
we could not observe any systematic tendency of [X/H] in terms of T c . 

These facts suggest that the "self-pollution" scenario proposed to explain the enhanced 
metallicity in planet-harboring stars, which invokes an accretion of metal-rich solid planetesi- 
mals (deficient in volatile elements while rich in refractory ones) onto the host star, is rather 
unlikely, since this mechanism would result in [X/Fe] < for [Fe/H] > for volatile elements 
and a systematic increase of [X/H] with T c . 

Therefore, as far as those 14 bright sample stars studied in this investigation are 
concerned, we would consider that the tendency of metal-richness in planet-bearing stars is of 
primordial origin (i.e., higher probability of planet formation for intrinsically metal- rich gas) 
rather than attributing it to an acquired character. 

This research have made use of the simbad database, operated at CDS, Strasburg, France. 
This work is based on observations carried out within the framework of the OAO key project 
"Comprehensive Spectroscopic Study of Stars with Planets", which aims to understand the 
property of planet-harboring stars by analyzing stellar abundances and line-profles, as well as 
to improve the precision of measuring stellar radial velocities to a level of planet-detectability. 
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Table 1. Data and adopted atmospheric parameters of the observed stars. 
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Note. The V magnitudes and the spectral types are from the SIMBAD database. The data of the mass, the half-major axis, 
and the perid for the planets around planet-harboring stars were taken from the "extrasolar planets" home page (http://cfa- 
www.harvard.edu/planets/encycl.html). The references for the T c ff , logg, [Fe/H] and £ values adopted in this paper are abbreviated 
in the last column: GL — Gonzalez and Laws (2000), G98 — Gonzalez (1998), DS — Drake and Smith (1993), GLTR — Gonzalez 
et al. (2001), GV — Gonzalez and Vanture (1998), FG — Feltzing and Gonzalez (2001), AC — Arribas and Crivellari (1989), F98 
— Fuhrmann (1998), and SHKTT — Sadakane et al. (1999). 

* Chromosphcrically active star with spectroscopic anomalies (cf. Drake, Smith 1993) 

t The companion mass may actually not be that of a planet, because its orbital inclination angle (i) is suspected to be very small 
(i.e., face-on system) according to Gatewood et al. (2001). 

* The value from high-excitation Fe t lines was adopted (because of being less affected by a non-LTE effect). 
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Table 2. Adopted gf values of important lines.* 
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74 


10 


.74 


-1, 


52 


i 


6156 


.76 


10 


.74 


-0, 


.93 


i 


6156 


.78 


10 


.74 


-0, 


,73 


Fe i 


6157, 


73 


4 


,08 


-1, 


,26 


i 


6158 


.15 


10 


.74 


-1, 


,89 


i 


6158 


.17 


10 


.74 


-1, 


.03 


i 


6158 


.19 


10 


.73 


-0, 


44 


i 


6300 


,30 





,00 


-9 


.82 


Zn i 


6362 


.34 


5 


,80 


+0, 


.15 


i 


7771 


.94 


9 


.15 


+0 


.32 


i 


7774 


.17 


9 


.15 


+0 


.17 


i 


7775 


.39 


9 


.15 


-0, 


.05 


N I 


8683, 


.40 


10 


33 


-0, 


.05 


S i 


8693 


,93 


7 


.87 


-0, 


,51 


S i 


8694, 


.03 


7 


.87 


+0 


.08 


C i 


8727, 


.12 


1 


.26 


-8, 


.21 


Fe i 


8727, 


.13 


4 


.19 


-3 


.93 


Si i 


8728, 


,01 


6 


.18 


-0, 


.61 


Si i 


8728, 


,59 


6 


.18 


-1, 


.72 


Fe i 


8729, 


.15 


3 


42 


-2, 


.95 


Si i 


8729, 


,28 


6, 


,18 


-2, 


30 



* All data were taken from the compilation of Kurucz and Bell (1995). 
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Table 3. Data of measured equivalent widths (in mA). 



Star 


C5380 


S6052 


O6300 


Zn6362 


07771 


07774 


07775 


N8683 


S8693 


S8694 


r Boo 


56.8 


35.9 




29.6 


158.8 


138.2 


125.2 


36.2 


52.0 


85.1 


v And 


41.5 


25.1 


4.6 


24.5 


122.9 


115.8 


96.4 


19.6 


29.1 


55.1 


16 Cyg B 


23.9 


11.2 


6.5 


18.8 


70.3 


60.7 


54.4 


10.8 


11.8 


30.8 


£ Eri 


9.8 


5 . 5 


5.2 


9.3 














47 UMa 


25.8 


13.8 


6.0 


22.8 


82.8 


79.5 


60.6 


13.0 


13.8 


33.2 


HD 89744 


43.0 


28.5 


4.9 


26.9 


137.3 


131.0 


111.3 


22.9 


29.9 


63.7 


51 Peg 


29.5 


18.2 


7.2 


25.4 














55 Cnc 


19.5 


20.1 


9.9 


26.6 


45.2 


49.9 


31.1 


10.5 


9.0 


20.3 


70 Vir 


14.9 


10.4 


9.5 


19.3 


62.9 


54.9 


42.2 


4.(1 


13.3 


19.5 


p CrB 


15.9 


10.8 


5.4 


19.1 


72.8 


73.0 


56.4 


(i.l 


11.4 


27.8 


HD 217107 


29.4 


23.6 


8.7 


26.8 














HD 52265 


37.1 


25.5 


7.4 


26.1 


118.8 


111.3 


88.0 


16.6 


27.4 


54.3 


HD 38529 


35.5 


27.0 


12.5 


36.2 


85.5 


82.8 


63.7 


17.0 


30.5 


42.0 


14 Her 


23.9 


23.7 


8.4 


29.2 














t Cct 


6.7 


4.8 


5.0 


10.7 














a CMi 


51.6 


24.9 


3.6 


15.9 


170.4 


164.1 


136.5 


30.2 


28.7 


59.3 


16 Cyg A 


24.6 


15.3 


7.4 


24.4 


79.6 


66.5 


59.4 


8.0 


16.0 


32.5 


S3 Vir 


38.7 


22.8 


6.6 


26.3 














Sun 


19.9 


11.2 


5.2 


19.7 


72.4 


62.7 


47.9 


6.4 


10.6 


27.3 
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Table 4. Comparison of the abundances from different lines. 



Star 


[C]s380 


[C]s727 


[o]S„ 


[O]6300 


[0] 7 773 


A7773 


[Si]6087 


[Si] 8728 


[S]6052 


[S]s694 




PNflBT 


[Fc]^ 29 


r Boo 


+0.31 


+0.42 


+0 


09 




+0.19 


(-0 


38) 


+0 


26 


+0 


38 


+0 


28 


+0 


59 


+0.18 


+0.21 


+0.27 


v And 


+0.19 


+0.29 


+ 


07 


-0.08 


+0.14 


(-0 


27) 


+0 


18 


-I) 


24 


-I) 


16 


+ 


28 


+0.10 


+0.09 


+0.11 


16 Cyg B 


+0.13 


-0.02 


+ 


06 


+0.07 


+0.07 


(-0 


11) 


+ 


09 


-I) 


10 


+ 


02 


+0 


08 


+0.05 


+0.05 


+0.09 


6 Eri 


+0.14 








+0.00 








-0 


04 






+0 


17 






-0.10 






47 UMa 


+0.10 


+0.05 


+0 


05 


-0.01 


+0.11 


(-o 


15) 


+0 


01 


+0 


05 


+ 


04 


+0 


08 


-0.08 


-0.05 


-0.01 


HD 89744 


+0.13 


+0.34 


+0 


02 


+0.07 


+0.10 


(-0 


31) 


+0 


26 


+0 


27 


+0 


16 


+0 


21 


+0.21 


+0.15 


+0.35 


51 Peg 


+0.26 




+0 


20 


+0.20 








+0 


16 






+ 


26 






+0.19 


+0.21 




55 One 


+0.33 


-0.34 






+0.37 


+0.29 


(-0 


06) 


+0 


26 


+0 


50 


+ 


69 


+0 


33 


+0.27 




+0.43 


70 Vir 


-0.17 


-0.11 


-0 


07 


-0.01 


-0.02 


(-0 


13) 


-0 


01 


-0 


02 


-0 


04 


+0 


15 


-0.08 


-0.06 


-0.08 


p OB 


-0.19 


-0.25 


+ 


08 


-0.23 


-0.02 


(-0 


15) 


-0 


13 


-0 


17 


-0 


09 


-0 


05 


-0.30 


-0.29 


-0.24 


HD 217107 


+0.36 




+ 


27 


+0.32 








+0 


27 






+0 


52 






+0.28 


+0.26 




HD 52265 


+0.15 


+0.33 


+ 


09 


+0.26 


+0.14 


(-0 


24) 


+0 


24 


-I) 


43 


+0 


21 


+0 


28 


+0.19 


+0.26 


+0.27 


HD 38529 


+0.31 


+0.25 


+ 


2:5 


+0.30 


+0.19 


(-0 


19) 


+0 


35 


-I) 


37 


-I) 


42 


+ 


60 


+0.37 


+0.41 


+0.39 


14 Her 


+0.46 








+0.35 








+0 


36 






+0 


70 






+0.39 






t Cot 


-0.13 








-0.18 








-0.21 






+0.03 






-0.58 






a CMi 


+0.15 


+0.26 


+0.07 


-0.16 


+0.17 


(-0.43) 


-0.05 


+0.04 


+0.00 


+0.07 


-0.08 


-0.17 


+0.07 


16 Cyg A 


+0.06 


-0.38 


+0.04 


+0.09 


+0.06 


(-0.14) 


+0.09 


+0.12 


+0.10 


+0.18 


+0.07 


+0.12 


+0.10 


Vir 


+0.14 




+0.10 


+0.04 








+0.16 






+0.11 






+0.11 






Sun 


8.67 


8.41 


8.87 


8.97 


8.83 


(-0.11) 


7.46 


7.f 


S4 


7.36 


7.19 


7.54 


7.62 


7.64 



Note. The presented stellar abundances are the differential ones ([X] = log ex — l°g e ©) relative to the solar abundance given in the 
last row (log 6q ; expressed with the usual normalization of logeH — 12.00) . The values in the 7th column are the non-LTE corrections 
(A = log enlte — log clte) for the O I 7771—5 triplet (mean of the three lines) showing an appreciable non-LTE effect, while LTE was 
assumed for the other two oxygen lines (O I 6158 and [O i] 6300) for which this correction is negligibly small (cf. Takcda 2001). 
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Table 5. Abundance variations to changes in the model parameters. 







v And 




A t At 

70 Vir 






(6140 K, 4.12, 1.4 km s -1 ) 


(5500 K, 3.90, 1.0 km s _1 ) 






(Analysis 


of individual lines) 




C I 5380 


-0.07 


+0.09 


-0.02 


-0.09 +0.10 


+0.00 


N i 8683 


-0.11 


+0.10 


+0.01 


-0.13 +0.10 


+0.02 


U I ooOU 


+0.04 


+0.12 


+0.00 


+0.03 +0.13 


+0.00 


O I 7771 


-0.11 


+0.05 


-0.08 


-0.16 +0.08 


-0.05 


C\ t 777/1 


-0.12 


+0.05 


-0.08 


-0.16 +0.09 


-0.05 


O i 7775 


-0.11 


+0.07 


-0.07 


-0.15 +0.10 


-0.04 


S i 6052 


-0.06 


+0.08 


-0.01 


-0.09 +0.10 


+0.00 


C t Ci5(TJ 

o I ODyo 


-0.08 


+0.08 


-0.02 


-0.11 +0.10 


-0.01 


S I 8694 


-0.08 


+0.07 


-0.04 


-0.11 +0.09 


-0.02 


ZjII I OOOZ 


+0.03 


+0.04 


-0.03 


-0.03 +0.07 


-0.03 






(6080-6089 A fitting 


analysis) 




Si 


+0.06 


+0.01 


+0.00 


+0.03 +0.05 


+0.00 


Ti 


+0.23 


-0.11 


+0.00 


+0.22 -0.04 


-0.05 


V 


+0.11 


+0.00 


-0.01 


+0.16 -0.01 


-0.03 


Fe 


+0.05 


+0.06 


-0.04 


+0.08 +0.04 


-0.07 


v.. o 


+0.15 


-0.06 


+0.00 


+0.10 +0.02 


-0.02 


Xi 


+0.08 


+0.00 


-0.06 


+0.07 +0.03 


-0.09 




(6156.5-6158.5 A fittin 


g analysis) 







-0.12 


+0.11 


+0.00 


-0.14 +0.10 


-0.01 


Fe 


+0.09 


-0.01 


-0.11 


+0.10 +0.00 


-0.19 






(8726-8730.5 A fitting analysis) 




C 


+0.01 


+0.11 


+0.00 


-0.05 +0.15 


-0.01 


Si 


+0.04 


-0.03 


-0.07 


-0.01 +0.01 


-0.07 


Fc 


+0.09 


+0.01 


+0.02 


+0.11 +0.00 


-0.03 



* AT eff = +150 K. 
f Alogg = +0.3. 

* A£ = +0.5 km s _1 . 
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Table 6. Final results of the abundances. 





rci 


[N] 


[O] [S] 


fZnl 




[Fc] 


[Co] 


[Ni] 


[V] 


[Til 


Slope 


a* 




75 


12(1 


180 648 


660 


1311 


1336 


1351 


1354 


1450 


1549 


with without 






Volatile elements 






Refractory elements 




X 


N 


Tb^, 


+0 35 


+0.47 


+0.14 +0.44 


+0 29 


+0 26 


+0.18 


+0.42 


+0.24 


+0.10 


+0 22 


-8.45 


-4.58 


v And 


+0 22 


+0.24 


+0.07 +0.22 


+0 07 


+0 18 


+0.10 


+0.26 


+0.12 


+0.14 


+0 08 


-2.78 


-0.76 


-lu v-^y& i-> 


+0.08 


+0.31 


+0.07 +0.05 


— 0.02 


+0.09 


+0.05 


+0.12 


+0.10 


+0.07 


+0.00 


-5.19 


+0.74 


e Eri 


+0.14 




+0.00 +0.17 


— 0.15 


— 0.04 


-0.10 


+0.00 


— 0.13 


+0.08 


— 0.05 




—8.16 


47 UMa 


+0.08 


+0.31 


+0.06 +0.06 


+0.06 


+0.01 


-0.08 


-0.01 


-0.04 


-0.05 


— 0.11 


— 15.81 ■ 


-10.77 


HD 89744 


+0.20 


+0.19 


+0.06 +0.19 


+0.17 


+0.26 


+0.21 


+0.38 


+0.16 


+0.27 


+0.24 


+7.99 


+9.32 


51 Pot? 


+0.26 




+0.20 +0.26 


+0.20 


+0.16 


+0.19 


+0.30 


+0.29 


+0.12 


+0.12 




—3.83 


55 Cnc 


+0.11 


+0.71 


+0.32 +0.51 


+0.51 


+0.26 


+0.27 


+0.55 


+0.50 


+ 0.46 


+0. 18 


— 3.06 


+6.15 


70 Vir 


—0.15 


-0.17 


-0.04 +0.05 


— 0.08 


— 0.01 


-0.08 


+0.02 


-0.09 


-0.05 


— 0.10 


+3.93 


+ 1.50 


p CrB 


-0.21 


-0.05 


-0.02 -0.07 


— 0.13 


— 0.13 


-0.30 


-0.15 


-0.30 


-0.21 


— 0.18 


— 9.42 


—8.32 


HD 217107 


+0.36 




+0.29 +0.52 


+0.31 


+0.27 


+0.28 


+0.49 


+0.41 


+0.33 


+0.20 




— 2.69 


HD 52265 


+0.21 


+0.15 


+0.14 +0.25 


+0.17 


+0.24 


+0.19 


+0.39 


+0.25 


+0.26 


+0.28 


+7.85 


+7.28 


HD 38529 


+0.29 


+0.46 


+0.23 +0.51 


+0.39 


+0.35 


+0.37 


+0.58 


+0.45 


+0.48 


+0.31 


+5.79 


+9.12 


14 Her 


+0.46 




+0.35 +0.79 


+0.56 


+0.36 


+0.39 


+0.67 


+0.59 


+0.58 


+0.35 




+0.11 


t Cct 


-0.13 




-0.18 +0.03 


-0.21 


-0.21 


-0.58 


-0.33 


-0.51 


-0.43 


-0.45 




-25.44 


a CMi 


+0.19 +0.27 +0.06 +0.03 


-0.12 


-0.05 


-0.08 


+0.10 


-0.13 


-0.07 


-0.25 


-19.31 


-15.45 


16 Cyg A 


-0.09 +0.04 +0.06 +0.14 


+0.11 


+0.09 


+0.07 


+0.16 


+0.14 +0.12 +0.09 


+7.40 


+7.85 


Vir 


+0.14 




+0.08 +0.11 


+0.10 


+0.16 


+0.11 +0.22 +0.13 +0.09 +0.10 




+ 1.77 


Sun 


8.58 


8.13 


8.87 7.28 


4.51 


7.46 


7.54 


4.72 


6.35 


3.89 


5.01 







Note. The abundance results presented here are the differential values relative to the solar abundances (given in the last row). 
* The slope (a; expressed in unit of 10 -5 dex K — 1 ) of the linear-regression line, [X/H]— a ■ T c + b. 
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Fig. la 



T 




£ (FeI+)TiI Nil ~ 

c/5 



6080 6082 6084 6086 6088 

Wavelength (A) 



Fig. 1. Spectral fitting at the 6080-6089 A region to determine the abundances of Si, Ti. V, Fe. Co, and 
Ni. Open circles denote the observed data, while the solid curves represent the best-fit theoretical spectra. 
Each of the spectra, arranged in the decreasing order of T e s, are vertically offset by 0.2 with respect to 
the adjacent one. 

(a) Higher T cff stars with 5750 K < T cff < 6500 K. 



Fig. lb 




6080 6082 6084 6086 

Wavelength (A) 



6088 



Fig. 1. (b) Lower T eS stars with 5200 K < T eS < 5750 K. 
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Fig. 2 




8726 



8728 8730 
Wavelength (A) 



8732 



17 

Fig. 2. Spectra fitting at the 8726-8730.5 A region to determine the abundances of C, Si, and Fe. A 
verticai offset of 0.1 is applied to each spectra relative to the adjacent one. Otherwise, the same as in 
figure 1. 



Fig. 3 




6156 



6157 6158 
Wavelength (A) 



6159 
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Fig. 3. Spectra fitting at the 6156.5-6158.5 A region to determine the abundances of O and Fe. A vertical 
offset of 0.1 is applied to each spectra relative to the adjacent one. Otherwise, the same as in figure 1. 



Fig. 4 
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Fig. 4. Comparison of the [Fc/H] values derived in this study with those taken from various literatures 
(cf. table 1). Filled circles correspond to planet-harboring stars, while open circles to the reference stars. 





20 

Fig. 5. Comparison of the abundances of five refractory elements with the Fe abundance: (a) Si, (b) Ti, 
(c) V, (d) Co, and (e) Ni. The same meanings of the symbols as in figure 4. 
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Fig. 6. Comparison of the abundances of five volatile elements with the Fc abundance: (a) C, (b) N, (c) 
O, (d) S, and (e) Zn. The same meanings of the symbols as in figure 4. 



Fig. 7 
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Fig. 7. Comparison of [C/Si] with [Si/H]. The same meanings of the symbols as in figure 4. 



